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Abstract: I discuss event cuts in DIS which improve the similarity between the 
current Breit hemisphere and hemispheres in e^e~ annihilation. I also present a 
method to study the scale evolution in quark hemispheres using data from fixed 
energy e+e^ experiments. This method can benefit from the high statistics and 
fiavour tagging at LEPl, also at scales relevant for a comparison with HERA results. 



1 Introduction 



In deep inelastic scattering (DIS), the current hemisphere of the Breit frame is expected to 
be very similar to one hemisphere in an e+e~ experiment. This expectation relies on the 
fundamental assumption of quark fragmentation universality. However, a set of features in DIS 
introduce corrections to this assumption. In this letter I will give an abbreviated presentation 
of the investigations in adressing two main issues: 

QCD radiation can give rise to high-p^ emissions, without correspondance in an e~^e~ 
event. A high-p^ emission in DIS sometimes manifests itself with a completely empty current 
region 0. A sizeable rate of such events introduces uncertainties to the interpretation of 
the data. Furthermore, problems arise even if the current region is not empty, but merely 
depopulated. It is therefore of interest to find a way to exclude high-p^ events using another 
signal than an empty current Breit hemisphere. In sections |^ and |^ we discuss an approach 
based on jet reconstruction. 

The flavour composition in e~^e~ and DIS is not exactly the same. When excluding high-px 
events from the DIS analysis, the boson-gluon fusion channel for heavy quark production is 
supressed. This implies a lower heavy quark rate in the studied DIS sample, as compared to 
e~^e~ data at corresponding energies. A uds enriched e~^e~ data sample with high statistics is 
available from the Z° peak. In section ||we discuss a method to study properties of e~^e~ quark 
hemispheres at different scales, using data from the fixed energy Z" experiments. 

2 Jet Algorithms 

In order to exclude from the DIS sample high-p^ events without correspondance in e"'"e~, it 
is natural to use fc_L-type cluster algorithms with a jet resolution set to Q/2. This scale is in 
analogy to the kinematical constraint for gluon emissions in e~^e~ , which is p±g < E^< 

In the HERA experiments, particles which in the lab frame have a large pseudo-rapidity 
w.r.t. the proton direction are not detected. In our analysis, we have chosen to exclude all 
particles with pseudo-rapidity larger than 3.8, to take this into account. For clustering purposes, 
an initial cluster is introduced along the proton direction, carrying the missing longitudinal 
momentum. 

Historically, k±_ algorithms were first designed for e^e~ physics, and a set of different al- 
gorithms exist 0. In DIS, jet observables depend on structure functions and an algorithm 
where the jet properties factorizes into perturbatively calculable coefficients convoluted with 
the structure functions is in general preferred. A algorithm carried out in the Breit frame 
and designed to fulfill jet requirements in DIS is presented in [Q, but we have chosen not to 
use it in the present investigation for the following reasons: 

We do not intend to investigate jet cross sections or the specific jets found by the clus- 
ter algorithm, but merely to exclude high-p^ events. In the accepted sample, the analysis is 
performed on the current Breit hemisphere, which is defined independently of any jet recon- 
struction scheme. Thus the infrared properties of the used algorithm, and in particular the 
problems adressed by the Breit frame algorithm in [Q, are not essential for the present study. 



1 




Figure 1: The fraction of MC simulated Figure 2: The relative difference between 
events with empty current Breit hemispheres mean energy in the current Breit hemisphere 
is significantly reduced after a cut in jet-px. and Q/2 is reduced after a cut in jet-p±. 



It is in stressed that the appeahng benefits of the algorithm requires the jet resolution 
scale Et to satisfy Aq,^^ -C Ef <^Q^ . By setting Et = Q/2, we would in the present study use 
the Breit frame algorithm in a way for which it is not intended. 

For these reasons, we have chosen to use the simpler e^e~ algorithms. To estimate the 
reliability of the results obtained using jets, we have used three different algorithms, LuCLUS 0, 
Durham [H and Diclus 0, applied in the hadronic CMS. 



3 Current Breit Hemisphere Properties 

To examine if a cut in jet-p^ is suitable to isolate a sample where universality is expected 
to hold, we study results for the current Breit hemisphere in MC generated events at HERA 
energies. The results are compared to e"'"e~ MC results. 

We simulate the electro- weak interaction in DIS using the Lepto MC [0. In both the 



e^e and DIS simulations, we use the Colour Dipole Model P], implemented in Ariadne |10 



to describe the parton cascade. The Lund string fragmentation model |]TT|, implemented in 



Jetset is then used to describe hadronization. The chosen MC programs give a generally 
good description of data in DIS and e+e~. 

The kinematical region of the DIS simulation is constrained by the cut > 4GeV^. For 
each Q^, the whole range of possible x- values is considered. Ideally, the properties of the current 
Breit hemisphere would depend on the scale only, but QCD radiative corrections introduces 
an implicit x dependence. When high-p^ events are excluded, the implicit x dependence is 
strongly reduced. To show that radiative corrections to the properties of the current Breit 
hemisphere are reduced by the cuts, it will suffice to consider the results integrated over x. 



Effects of the p^-cut 

In Fig |I|, the fraction of events in DIS with an empty current Breit hemisphere is shown to be 
significantly reduced after a p±< Q/2 cut, for all jet finding algorithms. 
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Another interesting observable is the average energy in the current Breit hemisphere, which 
in general differs from a corresonding e+e" hemisphere. Consider two hemispheres in the rest 
frame of an e+e" annihilation event. In general, they have different masses and hence different 
energies. In other words, a high-p^ emission in one hemisphere reduces the energy of the other, 
due to energy-momentum conservation. 

As for an e+e~ hemisphere, the energy of the current Breit hemisphere at fixed depends 
on emissions in this hemisphere and in the nearby phase space of the opposite hemisphere. The 
latter is very large in DIS at low x, but after our suggested event cuts, it is reduced by the 
condition p±< Qjl- The corresponding kinematical constraint in e^e~ is however |p| < Q/2, 
which is more restrictive. Thus the region for high-p^ emissions which can reduce the energy of 
the considered hemisphere is larger in DIS than in e^e~. Without reaching for a quantitative 
prediction, we expect the mean energy and multiplicity of the current Breit hemisphere to 
satisfy 

(i^cBH) < |, (A^cbh) < ^Ar,,(Q2), (1) 

where Nf.f,{Q'^) is the average multiplicity in an e^e~ experiment with invariant mass Q^. 

If 2 {Eq^yi) IQ would be much smaller than 1, the similarity between the current Breit 
hemisphere and an unbiased e^e~ hemisphere is poor, and it is not likely that reliable conclu- 
sions from a comparison can be drawn. On the other hand, if 2 (ii^cBH) /Q ~ 1? the current 
Breit hemisphere sample may be closely related to an unbiased e^e~ hemisphere. In Fig the 
relative energy shift 2 (-E'cbh) /Q ~ 1 is shown to be sizeable for the unrestricted event sample, 
but significantly reduced after imposing a < <5/2 cut. 



Flavour Compositions 

A cut in suppresses the boson-gluon fusion channel for charm production in DIS and reduces 
the charm rate. In our MC simulations, charm is suppressed by roughly a factor 2 in a p_L < Q/2 
event sample. A problem for the interpretation of this result is the fact that the MC simulation 
program significantly underestimates the rate of charm events. In the MC this is about 10%, 



while data indicate a much larger value around 25% |jT2| . However, a similar suppression ought 
to be expected also for real data. If boson-gluon fusion is a more important source of charm in 
nature than assumed in the MC model, the suppression may be even stronger. It would then 
be appropriate to compare with uds enriched e^e~ data. How to obtain such data at different 
energies will be discussed in section |[ 

Monte Carlo results for average multiplicities are presented in the left plot of Fig |^. The 
results for the low-p^ DIS sample are essentially equivalent for the three considered jet algo- 
rithms, and only the DiCLUS result is shown. In this MC generated event sample, heavy quark 
rates are well below 10%, and it is therefore better compared to a MC e^e~ sample with only 
uds. The agreement between the low-p_L and the e^e~ uds sample is better than between the 
unrestricted event samples also shown in Fig ^ 



Energy Scale Corrections 

As discussed previously in this section, high-p^ emissions close in rapidity to the considered 
hemisphere implies a reduction of the average energy in the current Breit hemisphere, (i^cBH), 
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Figure 3: MC simulations of multiplicities in e^e~ and current Breit hemispheres (CBH). 
LEFT: e^e~ results with s = Q^, allowing for all flavours (solid line) and only uds (dashed 
line), compared to CBH multiplicities using all events (crosses) and a jet-p± < Q/2 sample 
(diamonds). Both the exclusion of high-p± DIS events and heavy flavour e~^e~ events improves 
the agreement between the multiplicities. RIGHT: For low scales, the agreement is further 
improved if comparing the CBH multiplicity with e'^e' data with squared mass s = A (Ecbu)'^ ■ 



to a value slightly smaller than Q/2. It is then natural to compare the current Breit hemisphere 
multiplicities with e~^e~ data at a squared mass s = 4(i?cBH)^- The agreement between e~^e~ 
and Breit frame multiplicities is then significantly improved, as shown in the right plot of Fig |^. 
Only the DURHAM result is presented, but the other results are very similar. 

Though the focus here is on radiative corrections, we can not firmly exclude the possibility 
that there are differences between e+e" and DIS also in the hadronization phase, which could 
imply an non-perturbative flow of energy and particles between the current region and the target 
region. Such an effect is masked by the radiative corrections, and an unexpected breaking of 
quark fragmentation universality in the hadronization phase could hardly be seen in the average 
multiplicity. However, once expected differences are corrected for, unexpected features which 
break quark fragmentation universality can still be searched for in less blunt observables, like 
strangeness rates [O, energy spectra and higher multiplicity moments. 

Comparing the left and right plot in Fig |^, we see that the high-px cut has a relatively 
small influence on the average multiplicity, as compared to the energy scale shift. However, 
in Fig III we note that the rate of empty current Breit hemisphere events, which are excluded 
from the analysis more by necessity than by theoretical understanding, are reduced with the 
high-p_|_ cut. This is also the case for the energy shift, as seen in Fig ^ Thus we flnd that a cut 
in jet-p_L is a powerful step towards an event sample where quark fragmentation universality is 
expected to hold. 



4 Scale Evolutions in Fixed Energy e+e Annihilation 

A large sample of uds enriched events are available from LEPl at i/i = 90GeV. This is not 
the case for other energies corresponding to the HERA kinematical range. In this section we 
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Figure 4: Multiplicities in MC simulated e~ events including all allowed flavours (dotted line) 
are higher than for light quarks events (dashed line). The multiplicity obtained by the thrust- 
cone algorithm in a uds sample at the Z° pole (solid line) are in better agreement with the uds 
result for moderate scales between 4 and 8 GeV. This corresponds to current Breit hemispheres 
in DIS with between 16 and 64 GeV^. For larger scales the thrust-cone results differ from 
the uds sample in a similar way as the full flavour results. 

discuss how the scale evolution of e~^e^ uds hemispheres can be examined, using data from a 
fixed energy e~^e~ experiment. 

In an e~^e~ experiment at squared mass s, we consider an artificial scale Q^ax < ^- Three 
jets are reconstructed with a kj_ cluster algorithm, and events where p±> Qmax/2 are excluded. 
We study the multiplicity for particles where the rapidity, measured in the thrust direction, 
satisfies 



These cones in the thrust direction correspond to e'^e~ hemispheres with squared invariant mass 
Q^s^yr, and their evolution with this scale can thus be studied in a fixed energy experiment. (For 
a more detailed discussion, see e.g. p^ .) 

This "thrust-cone" method is designed to be similar to the Breit frame analysis. Two 
differences are present, due to the absence of a t-channel probe in e~^e~ . The scale and the 
direction which defines the quark jets are determined by the probe in DIS. In e~^e~, the scale 
Q'^s.x is chosen freely and the rapidity is measured along the thrust direction. It is possible 
to use the thrust direction also in the DIS analysis, which could improve the similarity with 
the e~^e~ thrust-cone algorithm. A thorough investigation of the quantitative effects of this 
adjustment is however better performed with a full detector simulation, as they may depend 
on how the thrust direction is reconstructed in DIS events with a large fraction of the target 
region undetected. 

In Fig ^ the multiplicity evolution of an e^e~ uds sample is compared to results including 
all flavours and results obtained by the thrust-cone algorithm in a fixed energy uds sample. 
The thrust-cone multiplicities are plotted as a function of the average energy scales obtained 
in the cones. These scales are somewhat smaller than the chosen Qmax/S, and the deviations 
are of order 5 to 10%. 



y\ > -^H^/Qm^^)- 



(2) 
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The aim of the thrust-cone algorithm is to reproduce the scale evolution of a e"'"e~ uds 
sample better than a sample including all flavours. As seen in Fig ^, this is achieved for 
moderate Qmax, between 4 and 8 GeV. This corresponds to current Breit hemispheres in DIS 
with between 16 and 64 GeV^, which is an important range at the HERA experiments. 

In rare cases where an emitted gluon gives the hardest jet, the thrust will be directed along 
the gluon momentum. Thus the average multiplicity in the thrust-cone analysis will get some 
contribution from gluon jets, which results in a systematic overestimation of the multiplicity. 
For Qmax > 8 GeV, the effect of gluon "pollution" is of the same order as the effect of the 
heavy quarks. However, in a realistic experimental analysis, the thrust-cone results applied on 
uds events at the TP peak will benefit from very large statistics compared to experiments at 
lower energies, and also from the fact that the scale evolution can be studied over a large range 
with the same detector. To conclude this section, our investigation indicates that it would 
be interesting to compare DIS data not only to full e^e~ results at different s, but also with 
thrust-cone results in uds samples from the Z° experiments. 

5 Summary 

The assumption of quark fragmentation universality implies that the current Breit hemisphere 
in DIS is expected to be very similar to a hemisphere in e^e~ annihilation. However, the ex- 
perimental situations are different, and several corrections to universality are present. We have 
here proposed event cuts to improve the expected validity of quark fragmentation universality. 

We suggest to exclude from the DIS analysis events with high-p^ emissions, which have no 
correspondence in e"'"e~ events. The p_L-scale can be reconstructed using jet cluster algorithms 
with the resolution scale Q/2. For these purposes, it is suitable to use k± clustering algorithms 
originally designed for analyses of e^e~ annihilation events. Using Monte Carlo simulations, 
we investigate three different types of cluster schemes, and find that the agreement between 
e^e~ and Breit frame results improves after a cut in jet p±, independently of the specific choice 
of e^e~ k± algorithm. 

In the accepted low-p_L DIS sample, heavy quarks are suppressed. This motivates a com- 
parison with uds enriched e~^e~ data, which are available from the experiments at the Z° pole, 
but not at lower energies. We have here presented a method, the "thrust-cone algorithm", to 
study scale evolutions of e~^e~ quark hemispheres, using data from fixed energy experiments. 
With this algorithm, uds enriched data with high statistics from the LEPl experiments can be 
compared to results for the current Breit hemisphere, over a large range of energies. 
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